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A brief overview of the methods commonly used to determine or estimate the black hole mass 
in quiescent or active galaxies is presented and it is argued that the use of mass-scaling relations 
is both a reliable and the preferred method to apply to large samples of distant quasars. The 
method uses spectroscopic measurements of a broad emission-hne width and continuum lumi- 
nosity and currently has a statistical la uncertainty in the absolute mass values of about a factor 
of 4. Potentially, this accuracy can be improved in the future. When applied to large samples of 
distant quasars it is evident that the black hole masses are very large, of order 1 to 10 billion Mq, 
even at the highest redshifts of 4 to 6. The black holes must build up their mass very fast in the 
early universe. Yet they do not grow much larger than that: a maximum mass of ~10^'' Mq is 
also observed. Preliminary mass functions of active black holes are presented for several quasar 
samples, including the Sloan Digital Sky Survey. Finally, common concerns related to the appli- 
cation of the mass scaling relations, especially for high redshift quasars, are briefly discussed. 



1. Introduction: Mass-estimation methods for active galaxies and 
quasars 

The Hubble Space Telescope has played a key role in our ability to detect supermassive 
black holes in the centers of nearby galaxies and to determine their mass through its 
high angular resolution. This has started a new, exciting era in which we are now able 
to study these massive objects and how they affect their environment and its cosmic 
evolution, and thereby get a more complete picture of structure formation and evolution 
over the history of the universe. 

There are different ways in which to determine or estimate the black-hole mass in active 
nuclei, and they all differ somewhat from the methods applied to quiescent galaxies. 
Therefore, to place these methods in perspective and to explain why the method of 
"mass-scaling relations" is preferred for distant active nuclei, I will start by giving a brief 
overview of the methods used to determine or estimate black-hole masses in quiescent 
and active galaxies. 

In galaxies containing a quiescent or dormant central black hole, its mass is measured 
by determining the velocity dispersion of stars or gas close enough to the black hole 
that their dynamics are dominated by its gravity (i.e., within the black hole's "radius of 
influence" ) by means of the virial theorem. While this appears the most direct way to 
measure the black-hole mass, this method cannot be applied to active black holes due to 
the strong glare from the central source. The most robust mass-determination method 
for active nuclei depends on its source type. Active galaxies are classified as Type 1 
and Type 2, depending on whether the broad emission-line region can be viewed by us 
(Type 1) or not (Type 2). Type 1 sources tend to be oriented more face-on, and Type 2 
sources are typically more highly inclined, such that gas and dust in the host galaxy 
or in the outskirts of the central engine block our direct view of the broad-line region. 
For Type 2 sources, water mega-masers emitted by the molecular gas disk in the galaxy 
(e.g., Miyoshi et al. 1995) is a very reliable and accurate way to determine the black 
hole mass. For less-inclined Type 1 sources, the mega-maser method cannot be used. 
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The most robust way to determine their black-hole masses is through the emission-line 
reverberation mapping method that relies entirely on the variability nature of the central 
engine (e.g., Peterson 1993; see also Peterson, this volume), and thus not on high angular- 
resolution observations as are needed for quiescent galaxies. At the present time, the 
monitoring data that has been collected allow the so-called one-dimensional reverberation 
mapping in which the broad emission-line time delay relative to continuum variations can 
be determined. This delay is the light travel time of photons from the central continuum 
region and is, hence, the distance of the broad line gas from the continuum region. The 
black-hole mass is determined from the virial theorem by also measuring the velocity 
dispersion of this variable gas. Ultimately, we would like to do a two-dimensional analysis 
which also yields a velocity-delay map. With the latter, we would be able to determine 
the black- hole mass zero-point independently of other methods. The one-dimensional 
method does not allow this, so at present we are limited to estimating the absolute-mass 
zero-point by other means. This is explained next. 

For quiescent local galaxies, the masses of their black holes display a strong, well- 
established correlation with the velocity dispersion of stars in the galactic bulge, far 
beyond the gravitational reach of the black hole (e.g., Ferrarese & Merritt 2000; Gebhardt 
et al. 2000a), the so-called A^bh o"* relationship. Active black holes also display such a 
relationship (Gebhardt et al. 2000b; Ferrarese et al. 2001; Onken et al. 2004), but given 
the unknown geometry and structure of the broad emission-line region, it is not yet 
possible to independently determine the zero-point of the rcilationship. However, if one 
assumes that the active black holes in the local neighborhood have essentially built up 
their mass, we can assume that the two Mbh-c* relationships are, in fact, one and the 
same, and thereby empirically establish the absolute zero-point offset of the mass scale 
of active black holes (e.g., Onken et al. 2004). The assumption is fair, since most of 
the active nuclei in the reverberation-mapping sample have black holes accreting at low 
levels, between 0.1% to a few percent of the Eddington rate (Peterson et al. 2004). 

It is not always possible to apply the primary methods of reverberation mapping or 
mega-masers to determine the black hole mass. This is where "secondary or tertiary 
methods" become useful. They typically build on or approximate the primary methods. 
Table 1 gives an overview of these methods, which I name "Secondary" for convenience. 
For example, BL Lac-type objects are viewed at such small inclination angles that the 
continuum emission is very strongly Doppler amplified, and the source variability reflects 
details of the powerful radio jets we are looking right into, rather than that of the 
continuum and line-emitting regions. Also, BL Lacs typically have very weak or no broad 
emission lines, most likely due to the strong continuum boosting. For these types of 
sources, researchers have utilized the fact that they are hosted by elliptical galaxies. The 
physical parameters of elliptical galaxies (surface brightness Eg, effective radius, Ve, and 
stellar velocity dispersion cr*) correlate such that they form a "plane," the Fundamental 
Plane (Djorgovski & Davis 1987; Dressier et al. 1987). To determine the central black- 
hole mass, one can measure Eg and Ve so to infer cr* from the Fundamental Plane, and 
next use the Mbh o"* to infer the black-hole mass (e.g.. Woo & Urry 2002; Falomo et al. 
2003) . Unfortunately, owing in part to the finite thickness of the Fundamental Plane, the 
inference from this plane relationship itself can add as much as a factor of 4 uncertainty 
(Woo & Urry 2002) to the mass estimates. This method can also be applied when only 
(host galaxy) imaging data are available. Alternatively, one can determine the luminosity 
of the galaxy (bulge) and from the Mbh--^^ (Bulge) relationship (e.g., Magorrian et al. 
1998), a sister relationship to the A'/bh cr, relation, infer the black- hole mass (McLure 
& Dunlop 2002). The lowest uncertainties that have been reported are a factor 3 to 4. 
Unfortunately, this method is prone to larger typical uncertainties when applied to active 
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Table 1. Secondary mass-estimation methods 



galaxies, owing to the difficulties in separating the luminosity of the galaxy bulge from 
the strong nuclear source (e.g., Wandel 2002). With spectroscopic measurements of the 
galaxy bulge cth,, the black-hole mass can be inferred from the Mbh^c* relation with a 
significantly smaller uncertainty (Table; 1). However, as the Call triplet lines AA 8498, 
8542, 8662 shift into the water vapor lines for z S 0.068, it becomes increasingly difficult 
and the measurement uncertainties increase (e.g., Ferrarese et al. 2001; Onken et al. 
2004). Some authors have then resorted to determining cr* from a collection of stellar 
features, although weaker, at shorter wavelengths (e.g., Barth et al. 2005). 

Because reverberation mapping requires a lot of telescope time and resources, this 
method is difficult or impossible to apply to large samples of sources, especially for distant 
objects (e.g., Vestergaard 2004a; Kaspi 2001). It is also an issue that more luminous 
sources, owing to their larger size, vary on longer time scales and with smaller amplitudes, 
making reverberation-mapping analysis challenging. Even for the nearby PG quasars of 
order 10 years of monitoring data are required to constrain the black-hole mass well 
enough (e.g., Kaspi 2001). Instead, one method that has been applied to obtain mass 
estimates of large samples makes use of the fact that the width of the [O III] A5007 line 
approximates the velocity dispersion of the stars in the galaxy (Nelson 2000; Nelson & 
Whittle 1996). The idea is to use the FWHM([0 iii]) as a substitute for cr* and then the 
-^-^BH 0-* relationship to infer the central mass. However, the la statistical uncertainty 
in this method is a factor of 5 (Boroson 2003). 

One other secondary method seems to perform better overall: the method sometimes 
referred to as "mass-scaling relations." The method uses measiiremcnts of line width and 
continuum luminosity from a single spectrum of an active nucleus to estimate the black- 
hole mass. It is an approximation to the reverberation-mapping method, and therefore 
is also based on the virial theorem. The method rc;lies strongly on another result from 
reverberation mapping, namely the radius-luminosity relationship, based on the size of 
the line-emitting region (as measured for a particular emission line) which can be esti- 
mated from the continuum luminosity. The advantages of this method are that it can 
both be applied to nearby and distant Type 1 active nuclei and to large samples with 
relative ease. The methods relying on accurate measurements of the host galaxy proper- 
ties are limited to relatively low redshifts below z of about 0.5 to 1.0, since more distant 
host galaxies are difficult to characterize (e.g., Kukula et al. 2001). Mass-scaling relations 
presently provide an accuracy of the absolute-mass values of a factor of about 4. This is 
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a method for which we can potentially even further improve the mass estimates in the 
future. 

In this contribution, I will focus on mass-scaling relations as my preferred method to 
obtain mass estimates of black holes in large samples of distant quasars: (a) because 
of their easy application to a larger redshift range than the other methods, without 
going into the infrared regime; and (b) due to the lower uncertainties associated. Before 
I discuss this method, I will address the key component of this method — namely the 
broad-line region radius-luminosity relationship (Section [21). 

Section [3] is dedicated to mass scaling relationships and important related issues. The 
distribution of quasar black-hole masses over the history of the universe is described 
in Section IH while the first mass functions based on various large quasar samples are 
presented in Section [3 Potential issues with mass-scaling relations as commonly high- 
lighted in the literature are discussed in Section [6l before the summary and conclusions 
(Section [ZD. 

A flat cosmology with Hq = 70 km s^^/Mpc and fix — 0.7 is used throughout. 

2. The radius luminosity relationship 

Peterson (this volume) has already introduced the radius-luminosity relationship: the 
relation between the characteristic size of the gaseous region emitting a particular broad 
emission line and the ionizing continuum luminosity that excites this gas. Therefore I 
will not discuss these relationships and their history in detail, but only emphasize a few 
issues that are important for their applications to estimating the black-hole mass for 
distant active nuclei. 

2.1. The optical R-L relationship 

Peterson discussed the efforts to establish the physically relevant relationship between 
the size of the H/3-emitting region and the ionizing optical continuum luminosity without 
the contamination from stellar light in the host galaxy of the active nucleus (Bentz et al. 
2006a). Keep in mind the reason that the first R-L relations published had a steeper slope 
is mostly due to the contamination from the host galaxy. The degree of contamination, 
even from nearby luminous quasars thought to entirely dominate their host-galaxy light, 
is in fact larger (Bentz et al. 2006a) than originally assumed (e.g., Kaspi et al. 2005). 
Therefore, the only radius-luminosity relation for i?(H/?) and La (5100 A) that should 
be applied is that of Bentz et al. (2006a), for which the slope is firmly established to be 
1/2, because this reflects the intrinsic, physical relationship. 

2.2. The UV R-L relationships 

Kaspi et al. (2005) established that a similar R-L relationship exists for continuum 
luminosities in the UV and x-ray regions using reverberation results and archival data. 
In particular, i?(II/3) scales with the monochromatic continuum luminosity at 1450 A 
as i?(II/3) cx La (1450 A)°-^'^. Restframe UV luminosities are most useful when studying 
distant quasars, since we can then use ground-based observations in the observed optical 
region. However, in this case it is not the i?(II/?)-L(UV) relationship that is relevant, but 
rather that pertaining to the UV emission line (Mg il or C iv) with which it is to be used 
(Section [3]). Prior to 2007, no such relationship could be established for Civ, mostly due 
to the unavailability of a sufficient data base (see Peterson, this volume). The inherent 
assumption in the mass estimates based on C IV at the time was therefore that the 
i?(C iv)-L(UV) relationship that had to exist, would be similar to that of the i?(H/3)- 
L(UV) relation (e.g., Vestergaard & Peterson 2006). The recent successes of Peterson 
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et al. (2005) to measure a C iv time delay (or lag) in the low-lmninosity Seyfert galaxy 
NGC 4395 and of Kaspi et al. (2007) to measure a C iv lag for a luminous quasar has 
finally extended the luminosity range (to over 7 orders of magnitude in luminosity) of 
previously measured (reliable) C iv lags, allowing us to confirm that the slope of the 
i?(C iv)-L(UV) relation (of 0.53) is entirely consistent with that for H/3. This confirms 
the previous assumptions made and shows that, as it happens, the recently updated 
mass-scaling relationships based on Civ (Vestergaard & Peterson 2006), which assumed 
i?(Civ) cx L(UV)°'^'^, still hold without any modifications needed. 

It is worth noting that an R-L relationship has not been established for the Mgil 
emission line, since only one reliable determination of the size of the Mg ii-emitting regions 
exists at present (Metzroth et al. 2006). Mgii is the one emission line that has typically 
not been targeted in monitoring campaigns. Its isolation in the spectrum from the UV 
high-ionization lines and the optical Balmer lines has significantly contributed to this fact. 
The only reliable Mgii lag is based on improved reverberation data of NGC 4151, which 
show that Mg ii is consistent with being emitted from a distance twice that of C IV for two 
epochs; yet within the uncertainties, the two lines are also consistent with being emitted 
from similar distances (Metzroth et al. 2006). Apart from measurement uncertainties, the 
limitation that no contemporaneous measurements of the H/3 and Mgii lags currently 
exists. However, new H/3-monitoring data of NGC 4151 do suggest that Mgii and II/3 are 
emitted at similar distances from the continuum regions (Bentz et al. 2006b). 

2.3. Do the R-L relations apply to distant luminous quasars? 

There are a few reasons why we can comfortably apply these R-L relationships to quasars 
beyond our local neighborhood in which they are defined. First, quasar spectra look 
remarkably similar (to first order) at all redshifts. This has been well demonstrated by 
Dietrich et al. (2002), who show in their Figures 3 and 5 the strong similarity of the quasar 
broad-line spectra across the history of the universe when binning in either luminosity 
or redshift. While this complicates a thorough understanding of the broad- line region, 
it is a help in this case, because it shows that for all quasars the broad-line regions 
are essentially the same — whether you look at a quasar in your backyard or one at the 
furtherest reaches of the universe (e.g., Barth et al. 2003; Jiang et al. 2007). If the broad- 
line regions are similar, they will also obey the same radius-luminosity relationships. 
Second, a key point about the Kaspi et al. (2007) results is that the quasar for which the 
size of the Civ-emitting region could be determined is located at a redshift of 3! Hence, 
these results also indicate that high-redshift quasars do have similar broad-line regions 
that respond to luminosity changes as do nearby active nuclei. Third, when we apply 
the R-L relations to more distant quasars, we are, contrary what is sometimes stated, 
not extrapolating beyond the luminosity range for which these relations are defined. The 
R-L relation pertaining to the optical continuum luminosity La (5100 A) and i?(H/3) 
is defined over 5 orders of magnitude in luminosity: from 10*^ erg s"'^ to 10*^ erg 
(Bentz et al. 2006a; see Figure 3 in Peterson's contribution in this volume). The UV 
R-L relation, relevant for R{Civ) and L\{IS50 A) or L>,(1450 A) — these luminosities 
can be used interchangeably (Vestergaard & Peterson 2006) — is defined over 7 orders of 
magnitude in luminosity: from 10^^-^ erg s~^ to lO''^ erg s~^ (Kaspi et al. 2007). Next, I 
will compare these characteristic luminosity ranges to the luminosities of known quasars 
of various surveys. 

Figure [1] shows the luminosity distributions of a few quasar samples, including the 
Bright Quasar Survey (BQS) at 2 ^ 0.5, the Large Bright Quasar Survey (shown in 
two redshift ranges: z ^ 0.5 and 1.1 ^ z ^ 2.9 that are discussed in later Sections), the 
color-selected sample from the Sloan Digital Sky Survey Fall Equatorial Stripe (Fan et al. 
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Figure 1. Distribution of bolometric luminosities for various quasar samples (labeled). The 
dotted lines mark the maximum luminosity for which the optical (applied at z < 0.9) and UV 
R-L relationships are defined. Very few quasars, even at high z, have luminosities exceeding 
these, so the use of the R-L relations is reasonable. Left: The Bright Quasar Survey (BQS), 
the Large Bright Quasar Survey (LBQS), and the color-selected sample presented by Fan et al. 
(2001), as labeled. Right: The Sloan Digital Sky Survey Data Release 3 quasar catalog. Vertical 
dashed lines indicate the subsets where difi'erent emission lines (labeled) are used for the mass 
estimates (left to right: H/3, Mgii, Civ). 



2001), and the SDSS Data Release 3 quasar catalog (Schneider et al. 2005). These sam- 
ples are discussed further in Sections |4] and [5l I converted the monochromatic power-law 
continuum luminosities to bolometric luminosities, LboI, in order to show the luminosi- 
ties of the individual samples in the same diagrams. The continuum luminosities are 
determined through a decomposition of the quasar spectra into a power-law continuum 
component, an iron-emission component using the best and most complete iron tem- 
plates in the UV (Vestergaard & Wilkes 2001) and the optical (Veron et al. 2004), and 
a line-emission component. The LboI values are estimated by a simple scaling factor ap- 
plied to the monochromatic luminosities of 4.62, 4.65, 5.8, 9.74, and 10.5 (Vestergaard 
2004b, and based on Richards et al. 2006) for the La(1350 A), La(1450 A), La(3000 A), 
I/a(4400 A), and La(5100 A) luminosities, respectively. 

The luminosity distributions in Figure [1] show that, as expected, more distant quasars 
are also more luminous. The lower flux limits of the different surveys uncovering these 
quasar samples are evident by the sometimes sharp cutoff of the distribution toward lower 
luminosities. In each panel I have indicated the maximum luminosity, in LboI units, for 
which the R-L relationships are defined for the optical continuum luminosity, L\{5100 A), 
and for the UV continuum luminosities, Lx{1350 A) and La(1450 A). They are 10"*^ and 
]^q47.65 gj-g respectively, given the bolometric scaling factors introduced above. It is 
clear from Figure [1] that very few quasars, even the distant, most luminous ones, have 
luminosities beyond the luminosity ranges where the R-L relations are defined. This 
holds both at low and high redshift. Statements that the mass-scaling relationships may 
not be valid for distant luminous quasars because of a need to extrapolate far beyond 
the validity of the R-L relations for these sources are clearly not supported. 
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2.4. The power of the R-L relationships 

The R-L relationships have tremendous power — providing they are rehable, of course — 
as discussed here and by Peterson (this volume), we have good reasons to believe that 
they are. They allow us to determine black-hole mass estimates of thousands of quasars 
on a relatively short timescale from a single spectrum of each source. This is in sharp 
contrast to the large amount of time and resources required to determine the mass us- 
ing reverberation mapping. However, R-L relationships do not marginalize monitoring 
campaigns and the reverberation method. In fact, reverberation mapping is still very 
important, because it provides the time lags and sizes of the various line-emitting re- 
gions and the more directly and more robustly determined masses to which the scaling 
relations are anchored. 

3. Single-epoch mass estimates based on scaling relations 

Scaling relations are approximations to the mass estimates based on the reverbera- 
tion mapping method, as discussed by Peterson in this proceedings. Thus, similar to 
reverberation-based masses, the single-epoch mass estimates rely on the virial theorem, 
where Mbh = V^R/G, G is the gravitational constant and R is the distance from the 
black hole to the gas with velocity dispersion, V. It is thus assumed that the broad-line 
emitting gas that is used to probe the black-hole mass is gravitationally dominated by 
the black- hole potential. It has by now been firmly established to hold for most, if not all, 
active galaxy broad-line emitting gas (Peterson, this volume; Peterson & Wandel 1999, 
2000; Onken & Peterson 2002; Metzroth et al. 2006). 

Specifically, these (mass-) scaling relations use the broad-line width as a measure of 
proxy for the velocity dispersion of the broad emission line gas. The line width is param- 
eterized either as the line dispersion triine (the second moment of the line profile) or the 
fuU-width-at-half-maximum, FWHM. As discussed later and by Peterson et al. (2004) 
and Collin et al. (2006), the line dispersion is the preferred line- width measure for its ro- 
bustness. The continuum luminosity is used to estimate the distance to the line-emitting 
gas via the important R-L relationship introduced in Section [2l since, clearly, no time 
delays can be determined without extensive multi-epoch observations. In summary, the 
mass-scaling relations take the form of either: 

Mbh = CI X FWHM(line)^ x , 

or 

Mbh = C2 X af^,,^ x L^ , 

where CI and C2 are the mass zero-points for each of the equations and the index P is 
essentially 0.5 for continuum luminosities (Section [2]). The specifics of the most recent 
updates to these relationships for C iv and H/3 are presented by Vestergaard & Peterson 
(2006), to which the interested reader is referred. 

The possibility of using measurements from a single observation of an active nucleus 
to estimate the mass of its central black hole was first introduced by Wandel, Peterson, 
& Malkan (1999) in the optical wavelength region using the 11/3 emission line and the 
ionizing luminosity estimated from photoionization theory. Vestergaard (2002) expanded 
this into the UV-wavelength region using the CiV emission line by calibrating single- 
epoch measurements of FWHM(C iv) and L(1350 A) to black hole masses obtained 
with reverberation mapping data. This was followed by work using FWHM(C iv) and 
L(1450 A; Warner et al. 2003) and FWHM(Mgii) and L(3000 A; McLure & Jarvis 2002). 
In the latter study, Mgii is assumed to be emitted from the same distance as H/3, such 
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that the Mg ii FWHM is substituting the H/3 hue width. With the availabihty of UV- 
scahng relations, it became more easily possible to estimate black-hole masses of distant 
quasars residing at the highest observable redshifts, and for large samples thereof. 

3.1. An empirical determination of the f-factor 

In 2004, the reverberation-mapping database underwent a major revision (Peterson et al. 
2004), which resulted in improved determinations of black-hole mass and the radius 
of the variable H/?-emitting gas for the 36 active nuclei monitored at that time. As 
summarized in that paper (and Peterson, this volume), an /-factor enters the virial 
mass estimates owing to the unknown geometry and structure of the broad emission-line 
region. The appropriate value of / is largely unknown, but is expected to be of order unity. 
However, Onken et al. (2004) empirically determined the average value of the /-factor to 
be 5.5 for the reverberation sample of active nuclei. For those sources with reliable (and 
completely independent) determinations of both black-hole mass and host-galaxy bulge 
stellar-velocity dispersion, Onken et al. made the assumption that intrinsically the M-a 
relationships for active and quiescent galaxies are the same. The distribution of black- 
hole masses and stellar-velocity dispersion for the active galaxies show the same slope, 
but an offset in the zero-point, suggesting that the assumed /-factor at the time was not 
representative. Prior to 2004, we had to make assumptions about the value of / so to get 
an estimate of the black-hole mass. The most simple kincmatical structure, an isothermal 
gas distribution, with equal gas velocities in all directions, was often adopted, yielding 
an / value of 0.75. It is, in fact, a very important point that the first values of / adopted 
were entirely arbitrary, but were accepted due to the lack of a more qualified guess at this 
time. The work by Onken et al. (2004) has made all other previously assumed values of 
/ obsolete. It is worth emphasizing, since this is often misunderstood, that this empirical 
determination of / = 5.5 is valid for random inclinations of the broad emission-line region 
with respect to our line of sight. The average value of / was determined such that some 
sources will have the mass either overestimated or underestimated, but on average for 
the reverberation sample, the mass zero-point is consistent with the M~a relationship 
for quiescent galaxies. I will discuss in Section 16.11 that the effects of source inclination 
can be minimized simply by using the line dispersion to measure the line width, or by 
making an /-factor correction to the mass estimates based on the FWHM. 

It is important to note that if one seeks to determine the black-hole mass of a single 
object based on a single-epoch spectrum — as opposed to the masses for a sample of 
objects for which the mass-scaling relationships discussed in Section 13.21 can be used — 
then one should use / = 3.85 along with the line dispersion to determine the mass (Collin 
et al. 2006). Since a single-epoch spectrum resembles the mean spectrum much more than 
the RMS spectrum, and given that the mean and rms profiles are often different, the 
appropriate /-factor is also different. In other words, to recover the reverberation black- 
hole mass of one of the reverberation-mapped active galaxies using the line dispersion 
measured from a single spectrum, / — 3.85 is needed. But, if only the FWHM of the line 
is available, one should use Eq. 7 of Collin et al. (2006) to determine the appropriate 
/-factor, which depends on the FWHM value. 

3.2. Improved mass-scaling relations for Hf3 and Civ 

Since the Peterson et al. (2004) and Onken et al. (2004) papers significantly revised 
the reverberation data on which the Vestergaard (2002) scaling relations are anchored, 
a revision of the latter was necessary. Using a larger database than the 2002 calibra- 
tions, Vestergaard & Peterson (2006) present updated and improved scaling relation- 
ships for both the CiV and H/3 emission lines, calibrated to the revised reverberation- 



M. Vestergaard: Black-hole masses of distant quasars 



9 



mapping masses. Equations using either the FWHM or hne dispersion are presented 
for Civ. Note that since the scahng relationships are cahbrated directly to the rever- 
beration masses — deemed the most accurate mass determinations of the objects in the 
reverberation-mapping sample — the /-factor is absorbed in the zero-points of the mass 
estimates. The Vestergaard & Peterson study shows that, although the scatter of the 
single-epoch mass estimates relative to the reverberation masses has decreased somewhat 
owing to the improved data base used, there is still some scatter remaining, amounting 
to a la statistical uncertainty in the absolute values of the mass estimates of a factor 3.5 
to 5, depending on the emission- line and line- width measure. This includes the estimated 
uncertainty in the absolute zero-point of the reverberation masses of a factor of 2.9 or 
less, as estimated by Onken et al. (2004). This means that there is a 68% chance that 
the mass estimate of a single object is accurate to within a factor 3.5-5, but there is a 
95% chance it is accurate to within a factor of about 6-7. 

The remaining scatter is possibly due to several factors. For one, the intrinsic scatter 
in the R-L relationship dominates the uncertainties in the mass estimates. That scatter 
is likely related to the detailed geometric and kinematical structure of the broad-line 
region that we are unable to account for. In addition, active nuclei do vary in luminosity. 
The R-L and virial relationships indicate that as the nuclear luminosity changes, so 
will the size of the line-emitting region, and its location will dictate the gas velocity for 
a given black-hole mass. As a result, the R and V values should slide along the virial 
relationship, R cx V^^^"^, as the luminosity changes. This has been clearly demonstrated 
for the Seyfert 1 galaxy, NGC 5548 (Peterson & Wandel 1999, 2000) and for several broad 
emission lines. And the mass estimates should always be consistent. However, deviations 
from a tight relationship may be expected in a couple of situations. For example, if the 
broad-line region has not yet adjusted to the luminosity change such that the actual 
R and the measured V may not accurately reflect the measured L. With BLR sizes of 
hundreds of days for luminous quasars (e.g., Kaspi et al. 2000, 2007) this can perhaps 
never be avoided. Another possibility for deviation is when the density distribution of the 
gas in the broad-line region is not smooth or "typical" (i.e., based on the R-L relation, 
which yields an ensemble average). The different emission lines could be emitted from 
regions lying closer or farther from one another than what is "typical." According to 
the Locally Optimized Cloud model of the broad-line region (Baldwin et al. 1995), the 
predominant part of the flux in a particular emission line is emitted from the part of the 
broad-line gas that is most efficient, at a given time, of emitting this line. The different 
emission lines will thus be emitted from different subregions in the broad-line region 
determined by the detailed physical conditions, such as the gas density and the incident 
photon flux. And since the broad line region is a dynamic place, the locations of most 
efficient C iv-line emission relative to that of H/3 may change in time and as the changing 
luminosity allows us to "probe" different part of the broad-line region. In other words, 
there may always be some scatter appearing in the R-L relationships, and thus always 
some intrinsic uncertainties in the mass estimates that will limit how well we can establish 
the black-hole mass, especially for individual objects. However, the hope is that as we 
learn more about the broad-line region, we are able to make appropriate adjustments 
to these relationships, so the accuracy of the black hole mass estimates and of distant 
quasars will improve. 

One additional factor that can affect mass estimates is the relative source inclination 
of the broad-line region to our line of sight. I discuss this separately in Section [6T] 

Unfortunately, the Vestergaard & Peterson calibrations were published before Collin 
et al. (2006) demonstrated that the H/? FWHM line widths may, in fact, provide a biased 
estimate of the broad-line gas velocity dispersion, depending on the inclination of the 
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active nucleus with respect to our line of sight. So, while the Vestergaard & Peterson 
scaling relations are anchored in the reverberation masses for which inclination effects 
are reduced, the relations based on FWHM do not directly take into account the possible 
inclination effects that can affect the measured FWHM values. This issue, for both H/3- 
and C iv-scaling relations, will be addressed in a forthcoming paper. 

3.3. Mass-scaling relations for Mgll 

The existing relations involving the Mgii emission line presented by McLure & Jarvis 
(2002) can no longer bo used because: (a) these mass estimates turn out to be inconsistent 
with those based on H/3 or CiV by up to a factor of 5 (Dietrich & Hamann 2004), and 
(b) the relationship has not been updated to the new reverberation-mass scale. I have 
instead established new relationships using a high signal-to-noisc subset of SDSS quasars 
for which Mgil and one other broad emission line (H/3 or C iv) is also present. These new 
relations make no prior assumptions of where Mg ii is emitted. The relationships and the 
details of this analysis will be presented in a forthcoming paper (Vestergaard et al. 2007, 
in preparation). 

3.4. Relationship involving the H0 luminosity 

For sources where the luminosity at 5100 A is not representative of the nuclear lumi- 
nosity due to contamination from the host-galaxy light or the high-energy tail of the 
powerful radio emission in radio-loud sources that cannot easily be corrected, the H/?- 
line luminosity can be used to approximate the size of the H/3-emitting region (e.g., Wu 
et al. 2004). This is valid because the H/3- line flux is directly proportional to the ionizing 
luminosity, which is the determining factor for the size of the broad-line region. Notably, 
since the C iv-emission line is partially collisionally excited, this assumption is not valid 
for C IV. This is also emphasized by the Baldwin Effect (an inverse correlation between 
the C IV cqiiivalcnt width and the continuum luminosity) : the continuum and C iv-line 
luminosities do not scale with one another. However, there is also no real need for an 
alternate luminosity measure — since neither the host-galaxy light, the high-energy tail of 
the radio emission, or even the UV iron emission between Ly a and C iv make a signifi- 
cant contribution to the 1350 A or 1450 A continuum luminosity (e.g., Elvis et al. 1994). 
A recalibration of the relationship involving the H/3-line width and H/3-line luminosity is 
presented by Vestergaard & Peterson (2006). 

This relationship is not applied to the quasar samples discussed in this contribution. 

4. Black-hole masses in distant quasars 

4.1. Quasar samples 

In this Section I discuss the results of applying the mass-scaling relations to large samples 
of quasars: namely, the Bright Quasar Survey (BQS; Schmidt & Green 1983), the Large 
Bright Quasar Survey (LBQS; Hewett et al. 1995), the Sloan Digital Sky Survey (SDSS) 
Data Release 3 quasar catalog (Schneider ct al. 2005), the color-selected sample of Fan 
et al. (2001) based on the SDSS Fall Equatorial Stripe, the z ~ 2 sample of radio-loud 
and radio-quiet quasars studied by Vestergaard et al. (2000) and Vestergaard (2000, 2003, 
2004b), and several z w 4 samples from the literature. These samples and the data, with 
exception of the LBQS and the SDSS DR3 Quasar samples, are discussed in detail by 
Vestergaard (2004b). Details of the analysis of the LBQS and SDSS DR3 samples will be 
presented elsewhere (Vestergaard et al. 2007, in preparation). However, it is appropriate 
to point out that the spectral measurements are made after modeling and subtraction 
of the Fell and Feiii emission and, in the case of the DR3 sources below z = 0.5, the 
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stellar light from the host galaxy; the LBQS quasars were selected to be strictly point 
sources, which limits the contribution from the host galaxy. The Mgii and H/3 lines, in 
particular, are embedded in broad bands of iron emission, which can be very strong in 
some sources. If the iron emission is not (well) fitted and subtracted, it will bias the 
line-width measurements. In particular, for Mgii one will systematically obtain smaller 
line widths because half of the broad Mg il-emission-line flux is blended with the iron 
emission (Vestergaard & Wilkes 2001). To minimize measurement errors due to noise, 
bad pixels, and narrow absorption, the line-width measurements were made on model 
fits to the broad emission line that are representative of the broad profiles. 

The black-hole mass estimates are based on the updated calibrations of Vestergaard & 
Peterson (2006) for II/3 and C iv; this work also lists the most current mass determinations 
and estimates of the BQS sample. For Mgii, the new relations discussed in Section [3] were 
applied. I am using the scaling relations based on the FWHM line widths, because most 
of the data are not of sufficient quality to allow the line dispersion to be significantly 
measured. Since the line dispersion is an integral measure, it is less affected by (residual) 
narrow-line emission, but is much more affected by flux in the extreme wings of the 
line profiles, as well as the placement of the underlying continuum level. For these data 
that are not of very high signal-to-noise, contrary to reverberation-mapping data, the 
use of the FWHM allows us to probe a much larger subset of the quasar samples. As 
discussed in Section 16.11 corrections can be applied to the mass estimates based on 
FWHM measurements to reduce potential effects of source inclination. While this is the 
best approach, it has not been done in the current work (for the sake of consistency), 
because a correction only exists for the H/3 line and not yet for the C iv and Mg ii lines. 

4.2. Does radio-loudness depend on black-hole mass? 

The sample selection of the z « 2 quasar sample makes it an ideal dataset for studying 
the properties of quasars with respect to radio-loudness, since each individual radio-loud 
quasar is matched in redshift and luminosity to a radio-quiet quasar in the sample. The 
distribution of masses in Figure [2] shows that when luminosity and redshift differences 
are eliminated, there is no difference in black-hole mass for quasars of different radio- 
loudness, as the BQS otherwise suggests (Laor 2000); that result is probably due to 
selection effects. 

4.3. Distributions with redshift 

The distributions of the black-hole mass estimates for the different quasar samples dis- 
cussed in Section are shown in Figure [21 The distributions at the low-mass end are 
strongly affected by the flux limits of the surveys by which the quasars were found. At low 
redshift, say z < 0.5, we can probe less luminous sources than at higher redshifts. Since 
the BQS is a "bright" quasar sample in the nearby universe, it has black-hole masses 
spanning 3 orders of magnitudes. It contains some quite luminous sources like 3C 273, 
which has a black- hole mass of a few billion solar masses. At redshifts higher than 1, we 
can typically only probe down to a few times 10® Mq. The lack of data points below this 
mass is therefore an artifact of our sample selection. We do expect the quasar population 
to extend much below the current flux limit cutoff. 

It is characteristic that for the quasars we do detect at high redshift, say at z > 1, their 
black-hole masses are typically very large — of order a billion solar masses or more. This 
holds even beyond a redshift of 3 to 3.5, at which the quasar space density is declining, 
and for the most distant quasars known at 2; > 6. This shows that black holes must build 
up their mass very quickly, since, for example, at a redshift of 6.3 the universe was less 
than 900 Myrs old. The first star formation is thought to occur at redshifts between 
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Figure 2. Distribution of black- hole mass estimates for different quasar samples. Tfiese mass 
estimates are based on the Vestergaard & Peterson (2006) single-epoch scaling relations and are 
computed in a flat cosmology witfi Ho = 70 km s^^/Mpc and fix = 0.7. A mass of 10^" Mq is 
indicated by the dotted line. Upper left panel: The Bright Quasar Survey (BQS), the intermedi- 
ate redshift sample studied by Vestergaaxd et al. (2000), Vestergaard (2003, 2004b), and various 
samples of ^; ^ 4 samples from the literature studied by Vestergaard (2004b). The measurement 
uncertainty (not including the statistical accuracy of the mass estimates) is indicated at the 
bottom of the panel. Upper right panel: The BQS, The Large Bright Quasar Survey, and the 
Color-Selected Sample in the Fall Equatorial Stripe (Fan et al. 2001). Lower panel: The Sloan 
Digital Sky Survey Data Release 3 Quasar Catalog. Vertical dashed lines indicate the subsets 
where different emission lines (labeled) are used for the mass estimates (left to right: H/3, Mgll, 
Civ). The error bars in each region indicate the typical uncertainty in the mass estimates and 
includes the statistical uncertainty in the mass zeropoint of about a factor of 4. The decreased 
densities of objects at redshifts of about 2.8 and 3.5 is due to the increased difficulty in selecting 
quasars at these redshifts, as in color space the quasar tracks overlap with the stellar locus. 



6 and 9, based on the presence of heavy elements in the spectra of quasars at z « 5 
(Dietrich et al. 2003a,b) and w 6 (e.g., Earth et al. 2003; Maiolino et al. 2005; Jiang 
et al. 2007; Dwck et al. 2007) and chemical evolution models (e.g., Friaca & Tcrlcvich 
1998; Matteucci & Recchi 2001). While the black holes may form earlier, this does show 
that given the short time available, supermassive black holes must evolve rapidly. 
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Figure 3. Distribution of Eddington luminosity ratios, I/soi/^'Edd for the SDSS DR3 quasar 

catalog. 

It is noteworthy that there appears to be a maximum black-hole mass for quasars 
at 10 billion solar masses. This mass value is marked in the diagrams of Figure [2] for 
reference. There are very few black holes with masses above 10^° Mq, even for very 
large samples like the SDSS. While SDSS does have an upper flux limit to the survey, 
these limits do not affect the quasar selection at high redshift. The number of objects 
with mass estimates above 10^" Mq are statistically consistent with the uncertainties in 
the black-hole mass estimates. Therefore it is not unreasonable to conclude that there 
appears to be a limit to how massive actively accreting black holes can become, and that 
this limit lies at or near 10^° Mq. At this point, it is unclear whether this limit is due to 
the dark matter halo in which the black hole resides, whether it has to do with details of 
the accretion physics, or whether the black hole is simply exhausting its accessible fuel 
supply. Of course, it could be due to a combination of all three possibilities. 

The distributions indicate that there is a decline in black-hole mass from a redshift of 
about 1.5 to the present. The SDSS survey does have a bright flux limit that deselects 
the brightest (and hence the most luminous) objects in the nearby universe. That is the 
reason that a,t z < 0.5 no SDSS quasars have a black-hole mass above about 10^ Mq at 
z « 0, while the BQS has a number of objects above this mass. Nonetheless, even for 
the bright quasar samples of BQS and LBQS, one can see a (rapid) decline in black-hole 
mass from a redshift of 0.5 to the present. This does not mean, of course, that black holes 
with mass above 10^ Mq do not exist in the local universe; the giant elliptical galaxy 
M87 is known to have a black-hole mass of a few billion solar masses. It simply shows 
that the activity of the most massive black holes rapidly declines from z — 0.5 to the 
present. They simply "turn off" and become increasingly quiescent. 

4.4. The Eddington luminosity ratio distribution with redshift 

With both the black-hole mass and the bolometric luminosities readily available, we 
can infer the Eddington luminosity ratio, I/Boi/^Edd, a crude estimate of the black- 
hole accretion rate. The distribution of Eddington ratios for the SDSS DR3 quasars is 
shown in Figure [3l as seen for the luminosity and mass distributions, the lower flux 
limits define the lower boundary to the distribution. As expected, most of the quasars 
appear to radiate between 10% and 100% of their Eddington luminosities. Even in the 
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Figure 4. Mass functions of actively accreting black holes at different redshifts and for different 
quasar samples, as labeled. The "SDSS" sample at a redshift of 3.6 to 5.0 is the color-selected 
sample of Fan et al. (2001). The SDSS DR3 mass functions axe preliminary. The downturn at the 
low mass end is not real, but caused by uncorrected incompleteness in these bins. The bottom 
panel shows how the mass functions shown in the top panels compare to each other. 

low rodshift {z < 0.5) universe there are many highly accreting quasars. Although some 
objects appear to have super-Eddington accretion rates up to a value of 3 (0.5 dex), 
this is not likely to be real. For one, the statistical uncertainties are entirely consistent 
with no significant population of sources with super-Eddington luminosities. In addition, 
the iBoi/^Edd estimate is based on the assumption of Bondi accretion (i.e., spherical 
accretion), which is not very likely to hold for a black hole fed by an accretion disk. 
Moreover, the ieoi/^Edd estimates are crude and have uncertainties at least of a factor 
of about 3.5 due to the uncertainties in the mass estimates. 

The distribution of the iBoi/^Edd values with redshift shows that at the highest red- 
shifts the quasars must be highly accreting to be detected. As a result, all known z > 4 
quasars are accreting at or near their Eddington limits. 

5. Mass functions of actively accreting black holes 

In the spirit of showing preliminary results from work in progress, I briefly present and 
discuss our first cut at the mass functions of actively accreting black holes of quasars at 
a range of redshifts. The mass functions for the different quasar samples are shown in 
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Figure ID Please note that the downturn of the mass functions at the low mass end is an 
artifact of incompleteness in these bins that is not yet corrected for. 

The upper left panel shows the BQS, LBQS, and the color-selected sample of Fan et al. 
(2001). The BQS (filled circles) and low redshift subset of the LBQS (filled triangles) 
show consistent shape and normalization. The high- 2: LBQS subset (open triangles) lies 
slightly higher than its low-z subset (at its high-mass end) and lies even higher than 
the color-selected sample (filled stars) located at higher redshifts yet. However, this is 
entirely consistent with the observed space density of quasars, which peaks at a redshift 
of about 2 to 3 and declines above and below these redshifts (e.g., Peterson 1997). This 
diagram also seems to indicate that the comoving volume density of 1 to 10 billion solar 
mass black holes at redshifts of 4 to 5 (the color-selected sample) approximately resembles 
that of the nearby bright quasars in the BQS (at z < 0.5). 

The upper right panel shows the mass functions for the SDSS DR3 quasar catalog 
divided into various redshift bins as labeled. For reference, the color-selected sample 
from the left panel is also shown. Again, we see a decrease in the normalization of the 
mass function beyond a redshift of 3. In the bottom panel, the individual mass functions 
from the top panels are overplotted to show their relative normalizations. In particular, 
the low-z DR3 subset (open diamonds) is entirely consistent with the other two low-z 
subsets of the BQS and the LBQS. A more detailed discussion of the finalized mass 
functions are presented by Vestergaard et al. (2009) and Vestergaard & Osmer (2009). 

5.1. Constraints on black-hole growth 

The ultimate goal is to use the black-hole mass estimates and the estimates of the mass 
functions of actively accreting black holes to constrain how black holes grow and evolve 
over cosmic time. This will be done by comparing these empirical results with theoretical 
models of how black holes can grow. While a more detailed analysis is yet to be done, 
it is interesting to do a simple comparison with one such theoretical model. Steed & 
Weinberg (2003) modeled how the distribution of black-hole masses change with redshift 
depending on the dominant growth mechanism. In Figure [5] I compare the distribution 
of quasar black-hole masses in SDSS DR3 with two such growth scenarios from Steed 
& Weinberg. The empirical data suggest that the black holes are more likely grow in a 
fashion in which the accretion rate depends on the black-hole mass rather than one in 
which the black holes obtain their mass in a short accretion phase, equal to the quasar 
lifetime. A more detailed discussion will be presented elsewhere. 

5.2. The space density of very massive black holes 
In Section|4]and Figure[2]it was shown that the typical black-hole mass at high redshift is 
a few billion solar masses and extends to 10 billion solar masses. Since the most massive 
black hole (in the elliptical galaxy M87) measured in the local neighborhood is only 
about 3 X 10^ Mq, and the measurement method does not exclude the determination 
of more massive black holes, one may question whether the single-epoch mass estimates, 
discussed here, may indeed be overestimates. I address this issue further in Section 
but it is instructive to also look at the high-redshift space density of black holes more 
massive than M87. Since the mass functions contain information on the comoving space 
density of the black holes as a function of their mass, I have integrated the current 
version of the quasar black-hole mass function above a mass of 3 x 10^ Mq for the 
subsets with quasars at and above a redshift of 4 (i.e., the mass functions shown by the 
open circles and filled stars in Figure [4|) . This measurement shows that such massive 
black holes are so rare that we need to probe the local neighborhood out to a distance of 
~290 Mpc from earth in order to detect just one such massive black hole. This is a volume 
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Figure 5. Distributions of masses of actively accreting black holes at different redshifts for 
the SDSS DR3 quasar catalog. The filled squares in each panel show the observed distribution 
with redshift for different luminosity ranges, as labeled. The filled squares show the mean value, 
while the vertical bars show the 10% to 90% ranges. The break luminosity, Lbrk, is the charac- 
teristic luminosity at the knee of the quasar luminosity function. Theoretical model predictions 
are indicated by open symbols: Open triangles: Short accretion phase build-up; Open circles: 
mass-dependent accretion build-up. 




about 25 times larger than what has been studied at present {R « 100 Mpc; Ferrarese 
& Merritt 2000; Tremaine et al. 2002). Owing to the uncertainty in the mass estimates, 
the mass functions may intrinsically be steeper than our estimates in Figure [4] show. In 
that case, the space density of the most massive black holes will be even lower, and an 
even larger local volume is needed. Strictly speaking, based on the uncertainties in the 
mass estimates, the mass functions could be shallower and thus the frequency of very 
massive black holes could be higher. We may be able to test this in the future. Figure 1 
of Ferrarese (2003) indicates that with a diffraction-limited 8 m telescope, in principle, it 
is technically feasible to detect 10^° Mq black holes out to distances of about 1500 Mpc. 
With a diffraction- limited 30 m we can possibly probe a volume of several thousand Mpc. 

In conclusion, it seems the explanation for the lack of very massive black holes in the 
local neighborhood is simply that of their extreme rarity and the relatively small volume 
that we have probed so far. However, as briefly discussed in Section 16. 4[ there are now 
suggestions that the brightest cluster galaxies may have black holes with masses of order 
10^° Mq. It is, however, still to be clarified whether or not these black holes are the 
quiescent equivalent of the very massive black holes found in distant quasars. 



6. Potential issues with scaling relationships 

There are a few issues that may potentially present problems for the single-epoch 
mass estimates, including those based on scaling relationships, in the sense that the 
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current estimates may in some sense be skewed or systematically biased. The issues most 
commonly mentioned in the literature include the effects of the relative inclination of 
the broad-line region to our line of sight, the potential danger of high- ionization outflows 
in quasars which may affect the C iv line emission, the unexpected narrow C IV profile 
for some active nuclei, and the question whether the single-epoch mass methods may 
overestimate the actual black-hole mass. 

6.1. Source inclination issues 

There are clear indications that the broad emission-line region may have the plane-like 
geometry because the measured gas velocities are larger for more inclined sources. This 
has been established for the H/? line by Wills & Browne (1986) and for the base of 
the Civ line by Vestergaard, Wilkes, & Barthel (2000) for radio-loud sources where the 
ratio of the radio-core flux to the radio-lobe flux yields a crude inclination measure: as 
the relative inclination of the radio jet increases, the core flux becomes less dominant 
relative to the radio lobe flux; notably, the direction of the jet is expected to be normal to 
the central engine's accretion disk. In both studies, the authors saw a lack of objects with 
small inclinations (i.e., near face-on) and large line widths. The largest line widths are 
only seen in the most inclined sources. This "Zone of Avoidance" distribution indicates 
that the gas velocities in a plane parallel or nearly parallel to that of the accretion disk 
are larger than perpendicular to this plane. 

Since active nuclei are expected to be randomly oriented (for Type 1 sources: within the 
inclination ranges that defines this subset: <45°; Barthel 1989), and the current mass de- 
terminations of active nuclei (both the reverberation-mapping method and mass-scaling 
relations) do not specifically account for source inclination, inclination effects can con- 
tribute to the uncertainties in the black-hole mass of individual active nuclei. For en- 
semble determinations, this effect is expected to average out owing to the way the mass 
zero-point was obtained (see Section l3.ip . 

Source inclination has been suggested to explain the scatter of the reverberation masses 
around the Mbh-ct* relationship (e.g., Wu & Han 2001; Zhang & Wu 2002). Collin et al. 
(2006) investigated this issue and found that this only appears to work in a statistical 
sense and is not supported by individual sources for which the source inclination is known 
or well constrained (see Peterson, this volume). However, Collin ct al. did establish that 
the FWHM measurement tends to be more sensitive to source inclination than the line 
dispersion. This is consistent with the finding of Peterson et al. (2004) that the line 
dispersion is the most robust line-width measure for mass determinations. Based on 
the reverberation-mapping data base, Collin et al. provide /-factor corrections to mass 
estimates based on the FWHM(H/3) values if the line-dispersion measurement cannot be 
obtained. This allows mass estimates for which the effects of the relative source inclination 
is decreased. 

However, this correction is currently only available for the H/3 emission line. Corrections 
applicable to the C iv and Mg ii emission lines will be addressed in a future publication. 

6.2. High-ionization outflows 

There are general concerns that active nuclei typically contain high-ionization outflows 
that potentially affect the C iv line emission, and therefore the mass estimates based 
thereon. These concerns are based on several observations to suggest this. Leighly (2000) 
finds for the subset classified as narrow-line Seyfert 1 galaxies (NLSls), that for increas- 
ing nuclear luminosity the C IV emission line displays an increasing blue asymmetry. 
Observations of NLSls in the far-UV (Yuan et al. 2004) support the interpretation that 
these sources have high-ionization outflows. IZw 1 is a typical example of a NLSl with a 
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broad blue asymmetric Civ profile. It has a FWHM(H/3) « 1200 km while its Civ 
hne is measm-ed to have FWHM of ~4400 km s^^ (Baskin & Laor 2005). Since NLSls 
typically have H/3 line widths below 2000 km s^^, Civ line widths several times that 
of the H/3 line are not very representative of the velocities of broad line gas dominated 
by black-hole gravity. This is the background for the discouragement of using C iv for 
mass estimates for NLSls, and in general if the Civ profile has a triangular and blue 
asymmetric shape similar to that of I Zw 1 (Vestergaard 2004b) . 

Because NLSls are luminous sources with highly accreting (albeit small) black holes, 
there is a concern that quasars that are also luminous and highly accreting may also dis- 
play such high-ionization outflows that would make mass estimates based on C iv biased 
or even invalid (e.g., Bachev et al. 2004; Shemmer et al. 2004; Baskin & Laor 2005). Based 
on the z > 1 samples shown in the upper left panel of Figure [2] I previously established 
(Vestergaard 2004b) that none of the high-z quasar spectra, at least in those samples, 
have C iv profiles resembling that of I Zw 1 . For the sources with the strongest asym- 
metries, I compared their luminosities and mass estimates with those of the remaining 
quasars and found no indication at all that the mass estimates were affected. The mass 
differences were at the 0.1 dex level and some subsets with asymmetries even showed a 
lower mass. Assuming that the quasar samples of that study are typical of the quasar 
population at those redshifts, it is fair to conclude that high-z quasars with Civ line 
profiles resembling that of I Zw 1 must be rare. 

However, it is still possible that the C iv emission line is associated with outflowing 
gas. It is well known that the C iv emission line is blueshifted relative to its restframe 
wavelength of 1549 A (e.g., Wilkes 1984; Espey et al. 1989; Tytler & Fan 1992). Based 
on about 3800 SDSS quasars, Richards et al. (2002) established that as the blueshift 
increases (up to a maximum of about 2000 km s"^), the line-equivalent width decreases, 
the profile shape changes, and the continuum luminosity increases. When comparing the 
line profiles of subsets of quasars binned by C iv blueshift, the authors concluded that the 
observed blueshift is generated by an increasing deficit of the red side of the profile. Their 
Figure 4 does suggest a significant broadening of the C iv profile with increasing blueshift. 
Nonetheless, to determine the effect on the single-epoch mass estimates, a comparison 
of the FWHM of the profiles with the most and the least blueshift is needed. Based on 
the entries in Table 2 of Richards et al. one can infer that, intriguingly, there is only a 
15% difference in the FWHM measurements, which translates into a mere 30% effect in 
the mass estimates. Regardless of the cause of the blueshift, since most quasars display a 
smaller blueshift than 2000 km s~^ and the current mass estimates based on C iv scaling 
relations have uncertainties of order of a factor of 4, the effects of the C iv blueshift is 
insignificant at the present time. However, as the mass estimates improve sufficiently this 
effect needs to be accounted for. 

6.3. The unexpected relative line widths of Civand Hf3 

Motivated by the blueshifts and potential blue profile asymmetries of C iv, Baskin & 
Laor (2005) investigated the C iv- and H/3-profile differences using non-contemporaneous 
optical- and UV-literature data of the BQS quasar sample. They find, in particular, that 
for objects with FWHM(H/3) > 4000 km s"\ the Civ profile is not broader than H/3. 
Existing monitoring data and the fact that C iv is a high-ionization line indicate that C IV 
is emitted from a region closer to the black hole than H/3. The Baskin & Laor study shows 
that FWHM(Civ) is not always larger than FWHM(H/?), and in fact the former rarely 
comply with being a factor broader, as expected based on the virial relationship, if 
H/3 is emitted from a distance twice that of C iv. As it turns out, this issue is much less 
severe when problem data such as low quality lUE data, strongly absorbed profiles, and 
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NLSls are eliminated from the database (see Vestergaard & Peterson 2006 for details). 
However, it does remain that FWHM(Civ) does not increase as fast as FWHM(H/3) 
and an inverse correlation between FWHM(C iv)/FWHM(H/3) and FWHM(H/3) is seen 
extending to ±0.2 dex in the FWHM ratio. Part of this difference can be attributed to 
measurement uncertainties. 

How well established is the factor 2 difference in the C iv and H/J lags? Of the few 
sources with reliable lags for both lines, three sources (NGC 3783, NGC 5548, and 
NGC 7469) display Civ lags that are about a factor 2 shorter than those of H/3 for 
the same epochs (Peterson et al. 2004). For NGC 4151 there are no contemporaneous 
epochs for the two lines, but the data are also consistent with this result (Bentz et al. 
2006b; Metzroth et al. 2006). But for two sources the measured Civ and H/3 lags argue 
otherwise. Existing data for Fairall 9 display a C iv lag that is 1.5 times larger than that 
of H/3 and has large errors (Peterson et al. 2004). Photoionization theory expects the 
Civ and Lya line gas to be emitted from similar distances, so judging from the Lya 
lag, a more reliable CiV lag may be at 2/3 the lag of H/3. For 3C 390.3 the Civ line 
is narrower than H/3 and the C IV lag is almost twice that of H/3but within the large 
errors, the two lines could be emitted from similar distances (Peterson et al. 2004). So 
judging by the most recent analysis of the reverberation data base of Peterson et al. 
(2004) there are good reasons to expect H/3 to be emitted from a distance twice that of 
Civ. But the data on Fairall 9 and 3C 390.3 suggest that perhaps not all sources have 
broad-line regions that are as neatly organized. As argued in Section [3. 2[ the broad-line 
region is a dynamic place and a neat onion-skin-like emission-line gas distribution is, in 
fact, not expected. This is supported by photoionization theory. In their Figure 3 Korista 
et al. (1997) present the strength of different emission lines as a function of gas density 
and incident ionizing luminosity. A comparison of the diagrams for C IV and H/3 shows 
a significant overlap of parameter space for which sufficiently strong H/3 and C IV can 
be emitted. Hence, their line regions are not necessarily mutually exclusive, and they 
can indeed be located closer than the factor 2 difference in distance. In fact, the non- 
contemporaneous nature of most UV and optical data being compared and the dynamic 
nature of the broad-line region, which also likely has a clumpy gas distribution, is likely 
part of the explanation for the C iv- and H/3-line widths not always scaling perfectly with 
one another. 

After all, the most compelling argument in favor of C iv-based mass estimates is that 
reverberation results show the virial products {RV'^) are consistent for all the measured 
emission lines, which includes CiV (Peterson & Wandel 1999, 2000; Onken & Peterson 
2002). 

6.4. Are the masses overestimated? 

When the scaling relationships are applied to samples of high-redshift quasars, the esti- 
mated typical black-hole mass is of order a billion solar masses and extends to a few tens 
of billions of solar masses (Section 2]). The validity of these large black- hole masses have 
been questioned (e.g., Netzer 2003) for a couple of reasons. Firstly, black- hole masses of 
10 billion solar masses or more are not found in the local universe among quiescent black 
holes assumed to have been active in the past. The most massive black hole detected in 
the local universe is that in the giant elliptical galaxy, M87, with a mass of 3 x 10^ Mq 
(Harms et al. 1994; Macchetto et al. 1997). Second, from the well-established Mbh-c* 
relationship one can infer that a 10 billion solar mass black hole will reside in a galaxy 
with a bulge that has stellar velocity dispersions above 400 km s~^ (e.g., Tremaine et al. 
2002). Such massive galaxies are not observed in the local universe. It is therefore fair 
to question whether the single-epoch mass estimates from the scaling relations system- 
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atically over-predicts the mass of the black hole by about an order of magnitude. This 
would bring the highest quasar black-hole mass estimates in full consistency with the 
black-hole masses in quiescent galaxies measured so far in the local universe. The meth- 
ods to determine their black-hole masses are by many considered more reliable since they 
rely on stellar velocities that are less prone to non-gravitational forces, such as radiation 
pressure, that is expected to be present and strong in the presumably violent environment 
near the accreting black hole. 

There are different considerations supporting that mass-scaling relationships do not 
systematically overestimate the black-hole mass. A full discussion can be found else- 
where (Vestergaard 2004b). Here I will summarize the main points. First, it is argued in 
Section 12.31 that the radius-luminosity relationship is fully applicable to high-z quasars. 
Hence, there is no indication that the use of this relationship should bias us in any way. 
Second, there is no indication that the kinematics of the broad-line region of high-redshift 
quasars are dominated by forces other than black hole gravity. The argument in favor is 
a combination of the fact that this is certainly found to hold for nearby active nuclei and 
quasars (Peterson & Wandel 1999, 2000; Onken & Peterson 2002; KoUatschny 2003) and 
quasar spectra look very similar at all redshifts (e.g., Dietrich et al. 2002). In addition, 
as noted in Section 12.31 the monitoring data of the z « 3 quasar for which Kaspi et al. 
(2007) was able to measure a C iv emission-line lag also support the notion that high-z 
quasar broad-line regions are very similar to those at lower redshift. Third, it is argued 
in Section 16.21 that high-ionization outflows are not important for a very large fraction 
of quasars. Fourth, one can demonstrate that it is unlikely that the mass estimates are 
systematically too large using results from photoionization modeling (e.g., Korista et al. 
1997). For the sake of this argument, assume an estimated black-hole mass of 10^ Mq 
that is a factor 10 overestimated. The virial theorem shows that a typical quasar with 
FWHM(Civ) of 4500 km s~-^ will emit the Ly a- and Civ-hnes from a distance of only 
^^33 light-days from the central source. At this location for a typical quasar luminosity 
of LboI of order 10**^ erg s~^, the ratio of photons-to-gas particles is so high that not 
only are these lines emitted inefficiently, but especially the Cm] -line emission cannot be 
generated. The prominence of the C ill] in quasar spectra argue for the presence of low- 
density gas subjected to lower ionizing flux. The key is that both the Lya and the Civ 
lines are much more efficiently emitted from these regions, that necessarily are located 
much further from the central source in accordance with the predicted distance from the 
R-L relationship. 

In conclusion, there are no obvious indications based on our existing knowledge of the 
broad-line region that the mass estimates based on scaling relations are systematically 
overestimated. Even the crudeness of using the FWHM of a single-epoch emission line 
that contains a contribution from non- varying emission line gas does not significantly or 
systematically affect the mass estimates (Vestergaard 2004b). 

In closing, it is interesting to note that recent findings show that a natural "satura- 
tion" of the bulge stellar- velocity dispersion occurs above cr, « 400 km s^^ (Lauer et al. 
2007). For the brightest cluster galaxies and other high cr, galaxies, cr, is thus not an 
accurate indicator of the central black-hole mass. On the other hand, the galaxy lumi- 
nosity indicates black- hole masses of order 10^" Mq. The apparent contradiction of the 
large quasar black-hole masses is therefore no longer that evident. Moreover, as I showed 
in Section [521 10^" active black holes are so rare that we should not expect to see 
any in the local volume that has been probed so far. 
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7. Summary and conclusions 

It is argued that mass-scaling relations are the preferred method to estimate black hole 
masses of distant active galaxies and quasars owing to the ease with which the method 
can be applied to large samples, its higher accuracy relative to alternative comparable 
methods, and the reliability of the mass estimates. Furthermore, the accuracy can po- 
tentially be improved. This is important because we seek to use these mass estimates 
to understand how black holes grow and affect galaxies and their evolution — those that 
host the active black holes and those in their neighborhood. At present, the statistical la 
uncertainty in the absolute mass values amounts to a factor of 3.5 to 5, depending on the 
emission line and line- width measure applied; the C IV relations have lower uncertainties. 

Distributions of black hole masses of distant quasars at a range of epochs is presented. 
The estimated masses are very large, a billion Mq or more, even at redshifts of 4 to 6; 
the black holes can clearly grow and mature very quickly in the early universe. However, 
a maximum mass of 10^° Mq is also observed, which shows that black holes reach their 
ultimate growth limit then and thus must shut off their activity. Preliminary black- 
hole mass functions of various large quasar samples is also presented and discussed; a full 
analysis will be presented elsewhere (Vestergaard et al. 2008; Vestergaard & Osmer 2009; 
Vestergaard et al. 2009, in preparation). Furthermore, I outlined and briefly discussed 
some potential issues associated with applying mass-scaling relations to luminous high- 
redshift quasars and concluded that none of these jeopardize this method or significantly 
affects the mass estimates. Also, there are no indications that the mass estimates are 
systematically too high. 

I am grateful for financial support of this work by NASA through grants HST-AR- 
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